Gravitational wave recoil in Robinson- Trautman spacetimes 
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We consider the gravitational recoil due to non-reflection-symmetric gravitational wave emission 
in the context of axisymmetric Robinson- Trautman spacetimes. We show that regular initial data 
evolve generically into a final configuration corresponding to a Schwarzschild black-hole moving with 
constant speed. For the case of (reflection-)symmetric initial configurations, the mass of the remnant 
black-hole and the total energy radiated away are completely determined by the initial data, allowing 
us to obtain analytical expressions for some recent numerical results that have been appeared in 
the literature. Moreover, by using the Galerkin spectral method to analyze the non-linear regime 
of the Robinson- Trautman equations, we show that the recoil velocity can be estimated with good 
accuracy from some asymmetry measures (namely the first odd moments) of the initial data. The 
extension for the non-axisymmetric case and the implications of our results for realistic situations 
involving head-on collision of two black holes are also discussed. 
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I. INTRODUCTION 

The possibility that a body recoils while emitting grav- 
itational radiation has been known for decades This 
problem has been considered in the literature by means of 
many approximated and semi-analytical methods as, for 
instance, the particle approximation [2], post-Newtonian 
methods^, and the Close-Limit Approximation [3], lead- 
ing to typical recoil velocities of few hundreds of km/s 
for some realistic cases. Such conclusions, however, 
have changed drastically due to some recent advances 
in numerical relativity [5[. In particular, recent numeri- 
cal simulations of the merging process of binary black- 
holes indicate that asymmetrical gravitational wave emis- 
sion can indeed induce the merger remnant to recoil with 
velocities up to several thousands of km/s. The physi- 
cal nature and possible implications of such considerably 
higher gravitational recoil are now under intense inves- 
tigation (see, for instance, [7]). The calculation of the 
recoil velocity as a function of the black-holes initial con- 
ditions is a particularly important hard task. Since the 
full non-linear regime of Einstein equations is extremely 
intricate and costly to analyze, some approximated or 
"empirical" formulas relating the recoil velocity and the 
initial data have been proposed^]. 

The Robinson- Trautman (RT) spacetime ^ is perhaps 
the simplest solutions of General Relativity which can 
be interpreted as an isolated gravitational radiating sys- 
tem and, hence, it is certainly pertinent to the study 



'Electronic address; 
t Electronic address: 



romacedoQfma. if . usp . br | 



asaa® ime . unlcamp . br | 



of the gravitational recoil effect. However, despite the 
many strong mathematical results on the RT solutions 
available in the literature, only a few exact examples of 
RT spacetime are indeed known in explicit form (see, 
for references, It is known, nevertheless, that a 

regular initial data, corresponding typically to a com- 
pact body surrounded by gravitational waves, will evolve 
smoothly according to the RT equation into a final state 
corresponding to a remnant Schwarzschild black- hole [IH, 
which can be at rest or moving with constant speed. Our 
aim here is to go a step further in the characterization of 
such final evolution state as function of the initial con- 
ditions. Our results are motivated and checked by some 
numerical analysis. The Robinson- Trautman partial dif- 
ferential equation has been analyzed numerically in the 
recent literature being particularly suitable to be nu- 
merically solved by means of spectral method s [ip. [3 . [lH . 
We will follow Oliveira and Damiao Soares[ij, [la] and 
adopt the Galerkin method 16] for our analysis. How- 
ever, as we will show, we will implement it in a different 
way that will allow us to get simpler equations and a 
better accuracy. 

The present paper has four Sections and one Appendix. 
In the next Section, the main aspects of axisymmetric RT 
spacetimes are presented briefly. We show, in particular, 
how to read from the final state of the RT evolution the 
mass and speed of the remnant black-hole. It is also 
shown that, as expected, for (reflection-)symmetric ini- 
tial data there is no radiation recoil. In such a case, 
the mass of the remnant black-hole and the total energy 
radiated away are completely determined by the initial 
data, allowing us to establish analytical expressions for 
the results about the total radiated energy obtained nu- 
merically in [13] and [l^]- Section HI is devoted to the 
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study of generic axisymmetric initial data. We show that 
a typical RT evolution can lead to a gravitational recoil. 
A Galerkin projection method is used to calculate the 
final black-hole speed. We show also how the final re- 
coil velocity can be estimated with good accuracy from 
some asymmetry measures of the initial configuration, 
namely the first odd moments of the initial data. In the 
last Section, we discuss the physical interpretation of the 
typical initial data considered in this work, emphasiz- 
ing their relation with the problem of frontal collision of 
two black-holes. The extension of our results to the non- 
axisymmetric case is also commented in the last section. 
The Appendix presents a direct proof of a mathematical 
result used in Section II, namely that, for regular initial 
data, the final state of the RT evolution does correspond 
generically to a Schwarzschild black hole moving with 
constant speed. 



II. AXISYMMETRIC RT SPACETIME 

The standard form of the Robinson- Trautman (RT) 
metric in the usual spherical radiation coordinates 
(m, r, 9, (j)) reads I IQ] 

2 

ds^ = -(k- 2— - r(ln Q"^)^ du^ - 2dudr + ^d^^, 
\ r J 

(1) 

where Q = Q{u,9,(j)), mg is a constant mass parameter, 
and dSl^ and K stand for, respectively, the metric of the 
unit sphere and the gaussian curvature of the surface 
corresponding to r = 1 and u = uq constant, which is 
given by 



^^QMl + ^V^lnQ^ 



(2) 



with corresponding to the Laplacian on the unit 
sphere. Vacuum Einstein's equations for the metric ([1]) 
implies the Robinson- Trautman non-linear partial differ- 
ential eauation[To| 



d f 1 



(3) 



In this paper, we will focus on axisymmetric spacetimes 
and hence we will assume hereafter that Q — Q{u,9). By 
introducing x = cos 9 one has 

K = Q^ + Q^[{l-x')Q,]^il-x')Ql (4) 



and 



where 



(5) 



[il~x^)K,]^ = {1-xY{QQ....-QL) (6) 
- 8{x - x^)QQ^^^ - 4{1 - 3x^)QQ^^. 



Integrating ^ and assuming a regular gaussian curva- 
ture K one has 



d 

du 



dx 



''{u,x) 



0, 



(7) 



implying that the quantity Qq — J^Q ^dx is constant 
along the solutions of ([5]). Notice that, from ([T]), the reg- 
ularity of the surface u and r constants precludes us of 
having Q = 0. The regularity of the gaussian curvature 
K, on the other hand, requires < Q < oo. We normal- 
ize our data in order to have qq — 2, implying that the 
area of the surface corresponding to r and u constants is 
always 47rr^ along the u-evolution governed by 

Several classical results assure that, given a geometri- 
cally regular initial data Q{0,x), the solution of ^ ap- 
proaches asymptotically a stationary {Qu = 0) regime. 
The stationary solutions of ([5]) are such that 



(1 



A 



constant, 



leading to 



K ~ Aarctanhx + B, 



(8) 



(9) 



where B is another constant. Regularity of K on the 
interval [—1,1] requires necessarily A — Q. On the other 
hand, Eq. ([4]) implies that the regular Q solutions for 
which K is constant are such that Qxx = (see the 
Appendix for a direct proof). Therefore, the stationary 
solutions of ([5]) are always of the form Q = a-\-hx with a 
and h constants. Nevertheless, our choice of go = 2 yields 

— = 1. We choose in this work a parametrization 
such that a = cosh a and b = sinha. 

Given a normalized regular initial data Q{0,x), the 
asymptotic solution of (O will be always of the form 
Q{oo,x) = cosh a -I- xsinha. The final configuration is, 
hence, completely characterized by the sole parameter a. 
In order to unveil its physical role, let us consider the 
Bondi's mass function [17| 



dx 



(10) 



which has several desirable properties to define an "in- 
stantaneous" mass for the solutions of ©, see, for in- 
stance, [14]. In particular, we have that M{u) > mo for 
normalized initial data and, for u — > oo, it reduces to 



M(oo) 



2 (cosh a 



dx 



: sinh < 



= mo cosh a 
mo 
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where v — tanh a can be interpreted as the final velocity 
along the z axis of the remnant black-hole 17 1 . 

The Bondi's mass pUj) corresponds to the temporal 
component of the Bondi's four-momentum, which for 
generic (non-axisymmetric) RT solutions is given bvflTj 



^"^("^ = 4^ 



i{u,9,^) 



dS, 



(12) 
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where S"^ is the unit sphere spanned by the usual coordi- 
nates 9 and and with area element dS, and a — 0, 1,2,3, 
with = 1 and 77* being the radial three-vector directed 
to the point (9, (j)) on the unit sphere. For axisymmet- 
ric configurations, the non- vanishing components of the 
Bondi's four-momentum are Po{u) — M{u) and 



Too / 2; 



(13) 



which corresponds to the momentum carried by the so- 
lution along the z axis. For normalized initial data one 
has for u —^ 00 



mo 



Vi — f' 



= (l,0,0,-«). 



(14) 



reinforcing the interpretation of v as the final velocity of 
the remnant black-hole. 

Notice that, for symmetric (even) initial data (5(0, x), 
Eq. ^ implies that the solutions Q{u,x) of ^ are nec- 
essarily even for w > 0, establishing that there is no grav- 
itational recoil {v = 0) in this case. Such a behavior is, of 
course, in full agreement with the expectation that grav- 
itation recoil should be due to non-reflection-symmetric 
gravitational wave emission. Therefore, for even situa- 
tions, the constraint ([7]) determines completely the final 
evolution state. 



A. Radiated energy: reflection-symmetric case 

The fraction of the initial mass A/(0) radiated away 
along the w-evolution governed by (O can be calculated 
exactly for even configurations. Following p^ . we define 



A 



M{0) - M{oo) 
iWXO) ' 



(15) 



which clearly corresponds to the fraction of the initial 
mass lost due to gravitational wave emission. For even 
configurations, w = and we have simply 



A = 1 



dx 



-1 QHo,x) 



(16) 



It can be shown that < A < 1. As an explicit example 
of this exactly soluble case, let us consider the first even 
initial data considered in the papers [H, [lB| , namely the 
prolate spheroid corresponding to 

Q\0,x) = Ql{l-e^x^), (17) 

with < e < 1. The constraint qq = 2 implies that 



Ql 



lln(i±i 
2e V 1 - e 



leading finally to 



■In^ 



1 -e 



(18) 



(19) 



This is the exact analytical expression for the curves ob- 
tained in |.14] and ^15,] from numerical simulations. For 
sake of comparison with the results of [Sin, Fig. □ de- 
picts a semi- log plot of A as a function of y = 1 — e, follow- 
ing their conventions. A very good agreement is found. 
One can proceed in an analogous way for any other even 
(reflection-symmetric) configuration, we will return to 
this issue in the last Section. The exact expression for 
A is certainly valuable to the investigation of statistical 
properties of the non-linear gra vitational wave emission 
as those ones considered in [ij, [l^. For instance, it is 
clear from (1191) that the non-extensive distribution func- 
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FIG. 1: Tlie fraction A of the initial Bondi's mass lost due 
to gravitational wave emission for the initial configuration 
(|17p . as a function of y = 1 — e. The curve is in very good 
agreement with that one inferred from numerical results in 
[15I]. Notice, however, that the non-extensive distribution 
function proposed in [3, [IB] is merely an approximation for 
y « 1, see pU)! . 

tion proposed in [3, EBl is only an approximation valid 
for small e. In fact, we have 



1 39 

A = -(1 - yr + ^i^-y)' + - y)'): (20) 

for y « 1 (or e « 0). Notice that A — > 1 for e — > 1. 



III. GENERAL SOLUTIONS 

The evolution of generic initial data Q(0,x) is a greater 
challenge. Since the gravitational recoil is clearly related 
to the odd part of the function Q{u,x), one might con- 
sider in first place some asymmetry measures of the initial 
data. The simplest ones correspond to their first odd n 
moments 



qn{u) 



_i Q'^{u,x) 



dx. 



(21) 



which obey —go < In < 90- For the generic final evolu- 
tion state Q{oo,x) = cosh a -|- xsinha, we have 
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^ ^ V / fc=2oddi \ / \ J / 



valid for odd n. Fig. [2] shows the first final odd moments 
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FIG. 2: Final odd moments gn(oo), n = 1,3,5, as functions 
of the recoil velocity v = tanha, as given by (|22|l . Notice 
that, for a given < luj < 1, one has |gi| > IqsI > Igsl > ■ ■ ■ 

qn{oo) as functions of the recoil velocity v. As we will 
show, the relevance of the first odd moments (pij) rests on 
the fact that one can construct, as a linear combination of 
them, a second approximately conserved quantity along 
the solutions of the RT equation ([5]) in the framework of 
the Galerkin approximation. 

A. The Galerkin method 

We introduce now a Galerkin decomposition for 
Qiu,x) 

N 

Q{u,x) = Y.^,{u)P,{x), (23) 

where Pe{x) stands for the Legendre_j)olynomials. By 
using standard projection techniques tl6!|, Eq. ([5]) can be 
written as the system of ordinary differential equations 

2£+l 

be = -^^iQ' [(1 - ^^)K.]^, Pi), £ = 0, 1, . . . , iV, 

(24) 

where the inner product is given by (/, — J_^fgdx. 
From ([6]) and (|23|) . one can see that the functions in- 
volved in the inner product in the right-handed side of 
(|24p are simple polynomials in x. The integration can 
be performed exactly for arbitrary N (with the help of 
algebraic manipulation software as Maple, for instance), 
yielding 5**^ order polynomials on the mode functions b^. 
Notice that here, in contrast to the approach adopted 



in [l4l Il5l | , no transcendental function is involved in the 
Galerkin approximation. Now, the Cauchy problem for 
the RT equation corresponds basically in choosing the 
initial value of the mode functions bi{u) according to 

bi{Q)^-^{Q{0,x),Pi), (25) 

and then to solve the Initial Value Problem (IVP) given 
by 1211). 

Equation (|24)) has some useful properties that are in- 
dependent of N. For instance, their stationary solutions 
{bi = 0) have necessarily bi = for £ > 1 and arbi- 
trary (constants) bo and bi. Indeed, for any regular initial 
data, the systems evolves into the final state Q{oo, x) = 
6o(oo)Po(2;) + foi(oo)Pi(x), with 60(00)^-61(00)2 = 1 for 
the normalized expected. The recoil velocity will 

be given simply by t; = 5i(oo)/6o(oo). Another useful 
property is that for an even initial data, one has 6^ (w) = 
for odd i and, consequently, v = Q. The accuracy of the 
Galerkin decomposition is determined by the truncation 
order N in ([23]). It can be controlled effectively here by 
checking the conserved quantity go along the w-evolution. 
Typically, the expansion with N Legendre polynomials 
in (|23p is accurate provided that max|6Ar(u)| be small 
enough. 

Finally, we are able now to consider the evolution of 
generic initial data. The recoil velocity v can be calcu- 
lated by solving the IVP corresponding to the system 
of ordinary differential equations (|24p with initial condi- 
tions (|25|) . The recoil velocity determines completely the 
final state for normalized initial data, allowing the study 
of any other relevant quantity as, for instance, the frac- 
tion A of the initial mass radiated away as a function of 
the non-reflection-symmetric initial data. 

We have performed an exhaustive numerical analysis of 
the system (l24l) . The considered initial data include the 
following simple but representative family 

Q{Q,x)^Qo{l + ax + l3x^ +-ix^), (27) 

where the constant Qo is always chosen in order to ensure 
the normalization — 2. Some particular elements of 
this family are presented in Fig. [31 

Fig. [4| depicts a typical evolution for the modes bg{u) 
governed by ([M)) for a particular case of the family (|27p . 
The recoil velocity v can be read from the final state of 
the evolution for any initial data. We notice that, for the 
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FIG. 3; Polar plot of some typical non-reflection-symmetric initial data Q{0,x) of the family (I27|l . The initial condition a, b, 
and c correspond, respectively, to the parameters a = 1/2, /3 = 1, 7 = 0;a = /3 = 0, 7 = —2/3; and a = 0, f3 — 4, -y — 3. The 
dashed lines correspond to the associated gravitational radiation content (without scale), see Sect. IV. 



family of initial conditions (|27p. we always have beiO) — 
for ^ > 3 and, in this case, iV = 8 is sufficient to assure 
typically an accuracy (controlled by the constant qo (u) — 
2) of the Galerkin approximation up to 1%. Initial data 
with high QniO) typically require higher N in order to 
attain a given accuracy. The radiation content of the 
initial data can also give some clues about the minimal 
necessary value of N, see Section IV. 



B. Estimation of the recoil velocity 

Despite that the IVP associated to the equation ((24|) 
can be solved with quite modest computational resources, 
an analytical estimation of the recoil velocity v from the 
initial data would be certainly valuable. Since the fi- 
nal state of the RT evolution is completely characterized 
by the sole parameter v for normalized initial data, a 
second conserved quantity besides qo would suffice to de- 
termine completely the final state and, consequently, to 
determine the recoil velocity v. Unfortunately, the RT 
equation ([5|) does not seem to have any other conserved 
quantity rather than qo . On the other hand, its Galerkin 
approximation (j24p does indeed have a second conserved 
quantity. Such a new conserved quantity, however, will 
be only approximately constant along the solutions of the 
full RT equation. Nevertheless, the approximation will 
be as good as the Galerkin approximation is accurate. 
In order to construct an explicit expression for the new 
constant, we remind that ([5|) implies that the moments 
ni) obey the equation 



where ai{u) are quadratic functions of the modes bi(u). 
For odd n, the inner product in ([^5]) will select only the 
odd-£ terms in x in the summation leading to the 

following linear relation between (j„(u) and ai{u) 



3mo(?„(w) 



2N 
odd I 



ae{u) 



1 



(30) 



The right-handed side of pO|) has exactly N terms, im- 
plying, therefore, that one can have at most N linear in- 
dependent equations of the type (|28|) . The linear relation 
between q„ and ai{u) given by (j30p involves a Hilbert- 
type matrix T§| and, in particular, it is always possible 
to find -|- 1 rational numbers ai such that 




aiq2i-i{u) = 0. 



(31) 



The quantity between parenthesis is conserved along the 
solutions of ([M]) and, therefore, it corresponds to our 
second conserved quantity. One could also truncate the 
summation in (I30p in a given obtaining partial linear 
combination of the odd moments that are constant along 
the solutions of ((24| up to deviations proportional to max 
|a£_|_2(u)|. The first of such partial linear combinations 
are 



&moqn{u) = {x'^,[{l-x^)K^]J. 



(28) 



(^ = 


0) 


91, 


(^ = 


1) 


91 


(^- 


3) 


91 



5 

393 
14 



21 



(32) 
(33) 

(34) 



From 



and (1231). we see that 



[(1 



2N 

E 



ai{u)x^ 



The coefficients in the above expressions and the ai of 
(29) (|3T|) can be calculated in a straightforward way by us- 
ing, for instance. Gauss elimination in (j30p . However, 



6 



1.2 




FIG. 4: Evolution of the modes be{u) governed by (|24p for 
the case (a) of Fig. [S] A'^ = 8 was used, leading to an accu- 
racy (controlled by the constant qo = 2) of 10~*. The final 
evolution state has bo = 1.0197 and hi = 0.20017 and, con- 
sequently, the recoil velocity is w = 0.19628 and the radiated 
energy fraction A = 0.05420, calculated according to 
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-1 -0.5 0.5 1 

FIG. 5: Plot of « X gi(0) for some typical regular initial 
conditions. The dotted line is the curve predicted by (|35|l . 
The detail depicts the plot of qi{0) x qi(oo). The assumption 
of gi(c)o) = qi{0) is, typically, a good approximation when 
qi(u) is the dominant moment. 



qi(oo) w 'Zi(O) and estimate the final recoil velocity v as 

iA_i^,„i±£). (35, 

V \ V 1 - V J Q (0, x) 

In particular, v has the opposite sign of (?i(0), see Fig. 
[5l We emphasize, nevertheless, that (|35|) will be accurate 
solely in the cases where qi{u) is actually the dominant 
moment. 



IV. DISCUSSION 

The physical properties of the initial conditions corre- 
sponding to the family (P7)) can be investigated by con- 
sidering their radiation content, which is determined by 
the (l/r)-decaying part of the Riemann tensor and is 
proportional to the quantityfl^, [l^ 



our numerical calculations show that, for the typical ini- 
tial data considered here, the first odd moment qi domi- 
nates over the other ones, implying that the typical vari- 
ations (q'i(O) — gi(oo))/(7i(0) are rather small. We no- 
tice also that the typical initial data of the family P7|) 
considered here has 191(0)1 > |(73(0)| > |95(0)| > • • •, 
in agreement with the magnitude of the odd moments 
for the final state, see Fig. [2l This situation can fail 
for some very specific initial conditions. For instance, if 
one has |g3(0)| > |gi(0)| > |g5(0)| > qi will vary 
considerably along the solutions of ([M]) . but the com- 
bination given by (|34p will be approximately constant, 
and so on. Fig. [5] presents numerical evidences con- 
firming these results. For practical purposes, whenever 
|gi(0)| > 1^3(0)1 > 195(0)1 > • • •, one can assume that 



D{u,x) = -{l-x^)Q^du[^y (36) 

With the help of (O and ([6]), one can show that for poly- 
nomials Q{u,x) in x, the function D{u,x) will be also 
polynomial in x. Moreover, D(u, x) is an even (reflection- 
symmetric) function for even Q(u,x). The dashed lines 
in Fig. 3 are polar plots without scale of |il'(0,x)| corre- 
sponding to the radiation content of the associated ini- 
tial data. The asymmetry in the gravitational radiation 
emission responsible for the final recoil is clear. We notice 
that initial data with larger max|D(0,a;)| will typically 
require a larger value of the truncation order N to attain 
a given accuracy in the Galerkin approximation. For in- 
stance, case (c) of Fig. 3 requires a truncation order 
larger than cases (b) and (a) to keep the same accuracy. 

Some cases of the family (l?fl) are specially interest- 
ing since they are good approximations for the Brill- 



Lindquist initial data[20j 

Q(0,x)=Qof^=^ + -=^) (37) 

which can be interpreted as the final stage (after the 
horizon merging) of a frontal collision of two black 
holes [2li [2^ . with the parameters 11 > and < w < 1 
related, respectively, to the mass ratio and to the infalling 
relative velocity of the two black- holes. The constant Qq 
must be chosen in order to assure go — 2. We have 

M,w) = 2 + +3— In . 

(38) 



For = I (the equal masses case), the function (|37)) is 
reflection-symmetric and, in this case, the final state of 
the evolution is completely determined by the constraint 
qo = 2. For ^ 1, one can estimate the recoil velocity 
for this head-on collision approximation by using (j35p . 
For the the initial data ([57)1 we have 



arcsin w 



1 / _ 1^ 1^ ^ _ ^ ^^^^ "-"^ ^- (39) 



For small values of w, the condition (l39l) reduces to 



= w. 



(40) 



Fig. 6 shows the dependence of v with w for some values 
of /i as predicted by p9p and some numerical results. A 




FIG. 6: Dependence of the recoil velocity v with the infalling 
velocity w of the two black holes with different masses, as pre- 
dicted by (I39p . and some results from numerical calculations. 
Notice that v —v if /i — > 1/fi. 

very good agreement is found again. It is interesting to 
notice that 



lim^gi(0) = 2^^ 



1 



(41) 



for the initial data ([57]) . implying from ([55]) that there 
exists a maximum recoil velocity for this configuration 



lim lul = 



< 1. 



(42) 



In fact, Eq. ([55)1 implies that v < w for any (see Fig. 
6), a behavior already noticed in the numerical analysis 
of [22]. One can also calculate the fraction ([2S]) of the 
initial mass radiated away for this case 



A = l- 



{u, w) 



%/l - w2 h{u,w) 
where v is given by ([5^ and 



(43) 



h{u, w) 



2(1 + /i^) 8i|^ + ^l^) 15(^2-^^-*) 4(/x-hAi^) +40// + ^i*) 



(1 



„2^2 



)^ (1 



1 



2w; 



In 



i + u; 
1 — If 



(44) 
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The aspect of the curves (|43|) are similar to that one 
depicted in Fig. 1. In particular, for small w, one has 

^=l %l%^'^ ^' + Oi^% (45) 

compare with (|20|) . Due to (|42p . one has A ^ irrespec- 
tive of ^ for w ^ 1. 

We finish by commenting that the non-axisymmetric 
case Q = Q(u,9,(j)) can also be investigated by means 
of a Galerkin method. For such a case, the Galerkin 
decomposition ([^5]) is based on the spherical harmonics 

N e 

Q(u,0,0) =^ ^ (46) 

and a system of equations equivalent to can be ob- 
tained. In this case, the stationary regime corresponds 
also to the case for that bim = for £ > 1. The constant- 
K final state will have the form 

Q{oo, 0, 4>) ~ boo + bio cos 6 + a sin cos -|- csin^sint/i, 

(47) 

where 6§Q—(&^Q-fa^-|-c^) = 1 for the normalized case. The 
non-vanishing coefficients now can determine the modu- 
lus and the direction of the Bondi's four-momentum and, 
consequently, the recoil velocity of the remnant. These 
topics are now under investigation. 
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APPENDIX A 

One can check easily by a direct substitution that 
Q{x) = a + bx, with — b^ = K, is a regular solu- 
tion of (HI. Nevertheless, a stronger result holds in this 
case: all regular solutions of (jl]) with constant K are nec- 
essarily of this form. We are interested in the geometri- 
cally regular solutions (0 < Q{x) < oo and |QK(a;)| < oo 
for —1 < X < 1). The relevant phase space is three- 
dimensional and spanned by {x,Q, Qx)- Notice that any 
solution such that Q^x — must be constrained on the 
surface C of the phase space corresponding to the points 
such that 

L{x, Q, Qx) = - 2xQQx - (1 - x^)Ql - K = 0. (Al) 
One can show that, along any solution of ([4]), one has 

{l-x^)Q^ = 2{xQ + {l-x^)Qx)L, (A2) 

confirming that Z, is indeed an invariant surface of 
The linear equation (|A2p has the solution 

Lix,Qix),Qx{x))^A^^, (A3) 

1 ~ X'' 

where A is a constant, implying that any solution of (|^ 
such that L ^ cannot be regular in x = ±1 
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